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ABSTRACT 

We study the Dirac equation of a charged massless spinor on the general charged AdS 
black hole of conformal gravity. The equation can be solved exactly in terms of Heun's 
functions. We obtain the exact Green's function in the phase space {(jj,k). This allows us 
to obtain Fermi surfaces for both Fermi and non-Fermi liquids. Our analytic results provide 
a more elegant approach of studying some strongly interacting fermionic systems not only 
at zero temperature, but also at any finite temperature. At zero temperature, we analyse 
the motion of the poles in the complex u plane and obtain the leading order terms of the 
dispersion relation, expressed as the Laurent expansion of uj in terms of k. We illustrate 
new distinguishing features arising at the finite temperature. The Green's function with 
vanishing uj at finite temperature has a fascinating rich structure of spiked maxima in the 
plane of k and the fermion charge q. 
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Introduction: Charged black holes that are asymptotic to anti-de Sitter spacetimes (AdS) 
play important roles in applying the AdS/CFT correspondence to study some traditional 
condensed matter physics, including superconductors and Dirac-Fermi systems. This paper 
reports that the analytic Green's functions in the momentum space of some strongly in- 
teracting Dirac-Fermi systems can be obtained from charged AdS black holes of conformal 
gravity. 

It was proposed that such a black hole is dual to some strongly coupled fermionic 
system at finite charge density such as non- Fermi liquids [Ij. The procedure of deriving the 
Green's function was spelled out in [21 [3], where the Dirac equation of charged spinor in 
the extremal Reissner-Nordstrom black hole was studied. Imposing the appropriate horizon 
boundary condition on the solution, one can read off the Green's function in the momentum 
space (w, k) from its asymptotic behaviour. Unfortunately, there is no analytic solution for 
the Dirac equation, and the numerical approach was adopted in Nevertheless, the 

Fermi surfaces, which are defined as poles of the Green's function with vanishing u, were 
determined, and the Green's function for small u on the Fermi surface was shown to take 
the form [21S|: 



G(uj,k) = : ^ , (1) 



where vp is the Fermi velocity, hi, h2,^kp, t^kp are constants; they were determined numeri- 
cally in [21 [3]. The system describes non- Fermi liquids for u^p < ^. It turns out that indeed 
the case with Vkp < \ can arise from the charged AdS black hole as well as v^p > |. 

While physics can be adequately studied numerically, it is more elegant and satisfying 
if we can find an analytical result. Recently, two equally-charged AdS black hole in L> = 5 
gauged supergravity was considered in [Ij. In some suitable extremal limit, the Green's 
function with u = can be solved exactly and the Fermi surfaces kp can be determined 
analytically. A procedure [2, 3J can then be used to derive the Green's function for small u 
and the analytical expressions of the constants in ([T]) were obtained [4j. However, the black 
hole of this example suffers from having a curvature singularity that coincides with the 
horizon. There also lacks an analytic expression for general G(w, k) and furthermore there 
is no analytic result for non-extremal black holes. The black holes considered in [21 131 H] 
are arguably among the simplest ones in usual two-derivative gravities or supergravities. It 
is thus unlikely to find new black holes without convoluted matter in these theories whose 
Dirac equation becomes exactly solvable. 

Charged black hole in conformal gravity: Charged AdS black holes can also arise 
in conformal gravity. Conformal pure gravity in four dimensions is a four-derivative the- 
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ory constructed from the Weyl-squared term. Its minimum coupling to the Maxwell field 
preserves the conformal symmetry. The Lagrangian is given by 



e 



= laC^^'P'^C^..^ + \aF^ , (2) 



where F = dA. The Weyl-squared term is crucial in constructing critical gravity [5] and 
supergravity |6]. It was argued that Einstein gravity can emerge from conformal pure 
gravity in the IR region [7]. Furthermore, conformal gravity can be supersymmetrised in 
the off-shell formalism [8]. 

The most general spherically-symmetric black hole up to a conformal transformation 
was obtained in [9]. For our purpose of studying the Dirac- Fermi system on a plane, we 
construct the torus-symmetric solution, given by 

ds^ = -fdt^ + + r^(dx'^ + dy^) , A = q(— --)dt, 
J V ro r / 

HV3), w w . 

J = [r - ri){r - r2){r - ro) , 

r 

= (^^2)^ + (rori)^ + (ror2)^ - rir2ro(ri + r2 + ro) . (3) 

Without loss of generality, we let r = ro > be the horizon. It follows that if (ri,r2) 
are real, neither is bigger than rg. They can also form a pair of complex conjugates. The 
thermodynamics of the black hole is analysed in |10j . Note that the solution contains three 
non-trivial parameters, associated with the mass, charge and the massive spin-2 hair. We 
set A = —3 so that the asymptotic AdS has unit radius. The solution becomes extremal if 
we let r2 = ro, (or ri = ro, but not both.) The solution of ri = was obtained in by 
employing the pseudo-supersymmetry related to off-shell conformal supergravity. 
Charged fermion and Green's function: We are now in the position of studying 
the Dirac equation for a charged massless spinor: 

Y{d, + \u:fTab-iqA,)^ = ^. (4) 

We follow the procedure outlined in [21 [3]. For static ansatz, the contribution of the spin con- 
nection in the Dirac equation can be absorbed by the scaling of the field ^ = {—gg"')^"^. 
In the Fourier mode $ ~ ^-luit+ikx^, t^\iIi momentum lying in the x direction only, the 
Dirac equation reduces, with suitable choice of P-matrices, to two decoupled equations for 
'I' = {1^1,1^2) ■ The equation for -02 is the same as that for ipi up to /c -f-T- —k. The ipaS are 
themselves two-component spinors. Let V^i = {ui,U2)^ and u± = ui i: iu2, we have 

u'^ + Xi{r)u+ = \2{r)u- , u'_ + Ai(r)M_ = A2(r)u+ , (5) 

i(w + ga) ik 
Mr) = f , A2(r) = y=. (6) 



(Note that if we had considered the massive Dirac equation, we would have instead A2(r) = 
^ - ^.) It follows that 

""++ Pi +P2(?')^i+ = , -u" +pi(r)u'_ +p2(»^)^i- = 0, (7) 
Pi(0 = -^, P2(r) = |Ai|2-|A2p+Pi(r)Ai + A;. (8) 

A2 

There are two independent solutions for ^' and they can be organised as in-falling and 
outgoing modes on the horizon. The nature of the back hole requires keeping only the 
in-falling mode, which is associated with the retarded Green's function. To determine the 
Fermi surfaces, one only needs G{uj = 0,k). The horizon condition for the a; = solutions 
is more subtle. In extremal black holes, it turns out that one mode diverges on the horizon 
and it should be dropped whilst the other mode converges. However, for non-extremal back 
holes, both modes are convergent. There is no rule to impose the horizon condition if one 
has only solutions with vanishing u. 

Once the horizon condition is appropriately imposed, the Green's function can be read 
off from the asymptotic behaviour of the wave solution, given by 

G, = -G-,^ = lim ^ = -i lim ^+ ~ . (9) 

r-!>oo m r-s>oo u+ -|- U- 

Here, Gi and G2 are diagonal entries of the 2 by 2 matrix of the Green's function associ- 
ated with the massless spinor. In particular Gi corresponds to the right-handed spinorial 
operator and the G2 corresponds to the left, which we shall not consider in this paper. In 
this "standard" quantization, the function ui is treated as a source and U2 is the response. 
It is clear that ui and U2 are symmetric in the wave function. This leads to an alternative 
quantization in which the roles of ui and U2 are reversed. In this case, the Green's function 
of the right-handed spinor is given by —G2 instead. 

We find that the wave equations can be solved exactly for our general black hole ([S]). 
The solution can be expressed in terms of Heun's functions, which are too complicated for 
detailed discussions on the fine properties of the Green's function. In what follows, we shall 
give some detail analysis for the simplest case: ri = and r2 = tq. We shall then present 
the salient results in the increasing order of the complexity of the black hole and more 
detailed analysis are presented in [TO] . 

Case 1: ri = and r2 = Tq. We begin with the simplest extremal solution, which was 
obtained in [H]. The general solutions for u± are 
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where ci and C2 are two integration constants, and W^^uiz) denotes the Whittaker function. 
Further notations are specified as follows 



5 + ^^ = ^/^-9^ z = f"^'' . , (11) 

^ >' '^0 ro(r — ro) 



and K* is the complex conjugate of k. As r approaches the horizon ro, we have z — >• — oo, 
and the Whittaker function behaves as m ~ e~2^2:'^(l + 0{z^'^)) . Thus, the in-falling 



solution, selected by the black hole, is given by (10) with C2 = 0. The retarded Green's 
function can now be read off straightforwardly, given by 



Gi = = i— , 7 = ^ , (12) 



where U is the confluent hypergeometric function, defined by 

^' ^) = T^r^^^TT^i^i(«! ^) + ^^ff^z''\Fi{a -b + l;2-b;z). (13) 
1 (a — + Ij i (a) 

We can now investigate the Fermi surfaces. For a; = 0, we find 

Gi = -G- = -./P^. (14) 
^ k-qro 

Thus it becomes clear that kp = qrQ defines the only Fermi surface for the standard quan- 
tization and kp = —qvQ for the alternative quantization. For both cases, we have = 0. 
Such an extreme situation of non-Fermi liquids were also observed in extremal Reissner- 
Nordstrom black holes PIE]- Our general G{u}, k) allows us to study the behaviour in greater 
detail. Let us examine the behaviour of the Green's function near the Fermi surface, which 
occurs at kp = qrQ. If lo/tq approaches zero much faster than (k — kp), we find that 

G, ^'^^(^ - =^ . (15) 

l 2ik-kp) 



{k-kp) + qJ'^u: + 2{k-kp){-'-^)y 



Comparing to the small-a; formula ([T]), we find 



Vp^ r^O-q^J^-^^^, /ii ~ - y/2kp{k - kp) , /12 ~ - 2(/c - A;^) ■ (16) 

Analogous result can also be derived in the alternative quantisation [10]. It is instructive 
also to consider how the Green's functions behave for small to after we have literally fixed 
k = kp. We find 

G{u;, kp) ~ -i + 2g(log(-— ) - 27 - x{^q)) , (17) 

where 7 is the Euler number and x is the digamma function. The divergence is logarithmic 
in small uj. 
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Our Green's function (12) for general uj and k enables us to find poles and even the 
dispersion relations in the large parameter regions of lo and k. Let us consider ro = 1 
and q = 1. In Fig. [T| we give two contour plots of |G| in the plane of complex frequency. 
The left has k = 3 and the right has k = 7. We see that the number of poles and zeros, 
corresponding to the white and dark dots, increases with larger k. There is a branch cut 
along the negative imaginary axis, and the right and left "curves" disjoin at the branch cut. 
Note that the gap between the left and right curves disappears when q approaches zero. It 
is also of interest to note that in the oscillating region with imaginary ly, we find that the 
small area of the lower left plane with « 1 contains an infinite number of zeros and 
poles. 




i 


1 







Figure 1: The contour plots of |G| in the plane of complex oj. All the poles (white dots) 
and zeros (dark dots) occur in the lower half plane. Their number increases for larger k. A 
branch cut in the negative imaginary axis creates a shift to the curves. 

We study a number of such contour plots and obtain the poles for a variety of k. This 
allows us to draw the motion of poles in the complex uj plane in Fig. [2j It is significantly 
different from the small- (w/ro) analysis [3J, in which the artificial curve connecting the poles 
for constant /c is a straight line. 

The data associated with the solid lines in Fig. [2] enable us to construct the dispersion 
relations for these curves. We find that a very good fit is to use the Laurent series expansion 
of u in terms of k: 

UJ = 1- + ao + aik H , (18) 
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Figure 2: This plot depicts the motion of poles in the plane of complex oj. The solid lines 
are the locations of poles for continuous parameter k. The dotted and dashed lines describe 
the artificial curves connecting the discrete poles for a given k, as in Fig. [TJ The number 
of poles increases for larger k. The negative imaginary axis is the branch cut. Analogous 
figure for the extremal Reissner-Nordstrom black hole, but only with uj/r^ « 1, was obtain 
in [3]. 

where aj are complex numbers. For example, for the curves I, I' and II, we have 



I 

r 
II 



a_i = 0.68 + 0.811, ao = -1.65 - 0.73i, ai = 0.97 - 0.07i , 
a_i = 0.05 + 0.051, ao = 0.83 - 1.69i, ai = -0.95 - 0.181 , 
a_i = 2.65 + 5.361, oq = -3.63 - 2.641 , ai = 0.93 - 0.181 . (19) 



These dispersion relations (18) are valid for Re(a;) > for the curves I, II and III, etc, and 
for Re(a;) < for the primed curves. In Fig. [sj we show that the constructed dispersion 
relations fit very well with the the actual poles at least for k up to 10. The coefficients 
a_i,ao and ai were obtained numerically. It would be of great interest to obtain their 
analytic expressions in terms of rg and q. 




Figure 3: The fitting of the dispersion relations ( 18 ) (solid lines) and the actual poles (dots). 



The left figure is for curve I and the right is for curve II. In both cases, the horizontal axis 
is k; the upper line is Re(a;) and the lower is Im(a;). 
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Case 2: ri = and r2 < Tq. This black hole was also obtained in [TT] and it is the 
non-extremal generalisation of Case 1. We label tq = r_(_, r2 = r_ and hence Q = r^r^. 
We shall not present the derivation but simply the Green's function: 

. 1-21(7 + 20) a; , 

G((jJ,k)=i ; — - , il = — — , 1 =-f-(r+ — r^). 

^ ' ^ 1 + 21(7-2^7) ' 47rr' ^ttV + J 

2Fi[l - -iq, l + i' -iq; § - 210; 



^ 2Fi[---i<z,--i5;^-2iO;-£-](4.r) ' " Vq ^^"^ 

Note that T is the temperature of the black hole. Taking the subtle limit r_ — )• r+, we 
recover the results in Case 1. 

For a given black hole specified by the inner and outer horizons (r_,r-|_), the Fermi 
surfaces can be obtained by the poles of G(0, k) for any given q. Several new properties 
emerge in the non-extremal case. In the extremal limit, the general solution of n_ with 
a; = is given by 

Thus we see that has to be real; otherwise, the solution becomes oscillatory on the horizon, 
implying instability. In fact, there is usually a decoupling limit in an extremal black hole 
so that the metric becomes its near- horizon geometry AdS2 x S^. The parameter u then 
measures the conformal weight of the dual operator on the AdS2 boundary, and hence it 
must be real. Furthermore, for v > 0, which we can now choose without loss of generality, 
we must have C2 = for the solution to be well defined on the horizon. Thus the horizon 
boundary condition can be fixed for the wave solution even with u = 0, for which the 
concept of in-falling and outgoing becomes mute. The situation is very different for the 
non-extremal case. There is no decoupling limit to obtain the horizon geometry x T^. 
The wave solution with cj = is not oscillatory regardless whether v is complex or real. 
In fact, as we shall see later, there is no apparent combination of parameters that can be 
recognised as ly in the most general back hole background. Thus for the non-extremal case, 
there should be no physical requirement that be real, and in our case it can be pure 
imaginary as well. Another new feature is that both modes of the w = solution on the 
horizon are non-divergent and hence there is no way to impose the horizon condition on 
the oj = solution alone. This adds an additional level of complexity for studying the 
Fermi-Dirac system at finite temperature. 

Another important new feature in the non-extremal case is that the Green's function 
G(0, k) is not real whilst it is always real in the extremal case. This implies that when we 
adjust the parameter k, it is unlikely that both the real and imaginary parts of 1/G(0, k) 
vanish simultaneously. Thus the maxima of |G(0, k)\ are in general not literally divergent. 
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However, we do find many examples of spiked maxima that are of orders > 10 . For 
all practical purposes, they can be viewed as divergent and the corresponding /c's can be 
regarded as Fermi surfaces. In Fig. [4| we present the plots of |G(0,A;)| and we see that 
spiked maxima occur; some of them are of the order 10^ and higher. For example, at 
kp = 43.18346539964 we have |G| ~ lO^^ for r± = ±1 and q = 50. 




Figure 4: The non-extremal black hole is specified by r± = ±1. The left plot is |G| as a 
function of k with q = 50. The right is a 3D plot of |G| with < A; < 30 and < g < 50. 

As in the case one, it is instructive to plot \G{uj,k)\ in the plane of complex u, which 
we present in Fig. [sj The black hole parameters are r+ = 2 = 1 /r_ and the spinor charge 
is g = 1. The pattern of zeros and poles becomes more complicated in the non-extremal 
case. An important new feature is that there is no branch cut in the non-extremal case. 
As we change the values of k, we find that new patterns of zeros and poles emerge with 
some isolated pairs spread in the complex plane. The motion of the poles in the complex 
plane and also the dispersion relation require further investigation in a future publication. 
Interestingly, we find that as the black hole approaches the extremality, some group of zeros 
and poles condense on the negative imaginary axis to form the branch cut. 

It is of interest to understand how in the extremal limit the quantity v is restricted to 
real values. Let r_ = r_|_ — 5. We consider < S/r^^- « 1. We find that if we let u to be 
pure imaginary, for a given q, no matter how small, as 5 becomes smaller, there will always 
be a Fermi surface. In fact there is an oscillatory factor 6'^ for imaginary z/ in the Green's 
function, implying that as 5 becomes smaller and smaller, more and more Fermi surfaces 
emerge for k lies between and qro, for which u is pure imaginary. The maximum kp for 
the Fermi surfaces becomes qr^ in the extremal limit. As 6 approaches zero, the number 
of Fermi surfaces for k < qr^ becomes infinite, signalling instability. This is the origin of 
the instability observed earlier from the wave function of imaginary i/ that has an infinite 
number of oscillations near the horizon of the extremal black hole. 
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Figure 5: In the non-extremal case, the patterns of zeros and poles of \G\ in the complex u 
plane become more complicated. The black hole parameters are r+ = 2 = l/r_ and q = 1- 
The left contour plot is for k = 1/2 with imaginary i' and the right is for A; = 4 with real u. 
There is no branch cut in non-extremal case. 

Case 3: ri < r2 = Tq. This is the general extremal solution with two non-trivial 
parameters. The general G{uj, k) is of confluent Heun's functions. The G(0, k) is simpler, 
expressible in terms of a hypergeometric function: 



Z = ^ +;q\^Fi[v-\q,\ + v-iq-l + 2u- ^] . (22) 

^ Vro(ro-ri) / ro - n 

Note that in the extremal case, we must impose reality condition on v and the Green's 
function is manifestly real. For large enough q, multiple fermi surfaces can emerge. Further 
discussions are given in |10j . 

Case 4: the most general black hole. As shown in [10], the general Green's function 
can be expressed analytically in terms of Heun's functions, namely G{ijJ, k) = ij^, where 

cM(o, /3*; \ - 2if)o, 1 - 2iJ^o, § - 2iJ^o, \ + 2i(Q + f^i); 

^ m[a, b; 0,H- 2il7o, ^ - 2i(Q + Qi); 

(ro - r2)ri P\ iQ n, 
a=7 ^ — ) ^ — ' V = — } ^7 ill — 



(ro-ri)r2' (ro - ri)r2 ' ro(ri - r2) ' (n - ro)(ri - r2) ' 

^ _ wr2 ^ _ urro _ uj 

^ (r2 - ro)(r2 - ri) ' ° (ro - ri)(ro - r2) 4ttT ' 

Vr2(ro-ri) 2r2(ro - ri) / 
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C = Ifl _ n\|-2i(Q+ni) Q _ r2)l+2i{Q-Q2) (n^) 

We shall not analyse the Green's function here; it suffices to mention that spiked maxima 
of |G(0, in the orders beyond 10^ do arise in this general case. Details can be found 
in [10]. It is important to note that there is no parameter combination in the above that 
resembles u. 

Conclusions: In this paper we present the T^-symmetric charged AdS black holes in 
conformal gravity. We report that the Dirac equation of a charged massless spinor in this 
background can be solved exactly in terms of Heun's functions. We give detail analysis 
for the simplest extremal black hole background and present some salient results of the 
most general situation, for which more details can be found in [10\. Whilst the near Fermi- 
surface behaviour in the extremal case is expected to be similar to the previous results, our 
general G{uj, k) allows us to study properties beyond Fermi surfaces in analytic way. We 
study the motion of the poles in the plane of complex lo and obtain the leading orders of 
the dispersion relation, expressed as the Laurent expansion of u; in terms of k, with the 
coefficients determined numerically. 

Furthermore, we have exact G{uj, k) for non-extremal black holes as well and there is 
a rich structure of spiked maxima in |G(0,A;)|. This allows one to study some strongly 
interacting fermionic system at finite temperature. We illustrated many new distinguishing 
features at finite temperature. Conformal gravity with the Maxwell field exists only in 
four dimensions. Our results may provide a new way of studying some two-dimensional 
condensed matter system such as graphene in the strongly coupled region. 
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